Mem. S.A.It. Vol. 76, 



Modelling NIR molecular lines for Miras 

W. Nowotny^ S. Hofner^,T. Lebzelter\ B. Aringer^ and J. Hron^ 

' Institut fiir Astronomie der Universitat Wien, Vienna, Austria 

e-mail: nowotny@astro .univie . ac . at 
^ Departement of Astronomy and Space Physics, Uppsala University, Uppsala, Sweden 



Abstract. The atmospheric structure of Mira variables is considerably influenced by pulsa- 
tion. Molecular absorption lines in the near-infrared (NIR), especially second overtone CO 
lines, show therefore a characteristic behaviour in time-series of high-resolution spectra. 
We computed synthetic CO line profiles based on a new dynamic model atmosphere and 
derived radial velocities (RVs) from the Doppler shifted lines. For the first time, we could 
quantitatively reproduce observations of the very typical, discontinuous RV curves. 
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1. Probing atmospheric {cinematics 

Mira variables are luminous late-type gi- 
ants, representing stars of low- to interme- 
diate mass in the late evolutionary stage 
called Asymptotic Giant Branch (AGB). Due 
to the low effective temperature (<3500K), 
molecules can form efficiently in the outer 
layers and dominate their spectral appear- 
ance. Towards the end of their evolution, such 
stars become unstable against radial pulsa- 
tions. Eventually, this leads to the condensa- 
tion of dust grains in the cool atmospheres and 
stellar winds develop. 

High-resolution spectroscopy in the NIR 
has proven to be a valuable tool to inves- 
tigate atmospheric kinematics. RVs derived 
from the wavelength shifts of spectral lines 
provide clues on the gas velocities throughout 
the extended, dynamic atmosphere. In the past, 
CO lines of the A v=3 vibration-rotation bands 
at/}a!l.6//mhavebeenusedto study deep pho- 
tospheric layers in Miras, which are heavily af- 



fected by the pulsation. As the emerging shock 
front passes through the region of line forma- 
tion, it leads to a very characteristic behaviour 
of the CO line profiles (line doubling around 
phases of maximum luminosity) and results in 
discontinuous, S-shaped RV curves during the 
lightcycle (e.g. Hinkle et al. 1984). 

In Fig.^ (left), a compilation of observa- 
tional data (extracted from FTS spectra) is 
shown, suggesting a similar picture for all 
Miras, independent of their parameters (spec- 
tral type, period, metallicity, etc.). Realistic at- 
mospheric models should reproduce this. 

2. IVIodelling line profiles and 
synthetic radial velocities 

Dynamic model atmospheres are needed to re- 
produce the complex phenomena within the 
outer layers of Mira variables. Recently, we 
presented synthetic line profiles based on mod- 
els from Hofner et al. (2003). A special radia- 
tive transfer code was applied, which also takes 
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Observed radial velocities - all Miras, CO Av=3 Synthetic radial velocities - CO Av=3 




phase ij)^ phase ifi^^^^ + 0.3 

Fig. 1. Left: RV curves of second overtone CO lines for a sample of typical Miras (data taken from Lebzelter 
& Hinkle 2002). The grey, thick line was fit to guide the eye. Right: RV curves for the same lines as derived 
from model calculations (R=70.000). Model A was used in Nowotny et al. (2005a,b,c), while Model B is 
an improved version, which yields more realistic results, especially a comparable velocity amplitude of 
ARV=20.6km-s-'. 



velocity effects into account. Results can be 
found in Nowotny et al. (2005a,b,c). It could be 
shown, that state-of-the-art models can quali- 
tatively reproduce the behaviour (line profiles, 
RVs) of spectral lines originating at different 
depths. 

Still, some differences remained compared 
to observations. For CO Av=3 lines the line 
splitting appears too weak and the extracted ve- 
locity amplitude ARV too low (cf. Sect. 6.2. in 
Nowotny et al. 2005c). Here, calculations are 
presented using a new dynamic atmospheric 
model, for which the parameters' were some- 
what tuned to fit this aspect. Line profiles for 
the CO 5-2 P30 line were computed for vari- 
ous phases during a pulsation period. The spec- 
tral synthesis was carried out exactly as de- 
scribed in detail in Nowotny et al. (2005b, c). 
As also found there, the RV curves look more 
complicated for the highest resolution calcu- 
lated (300.000), but become similar to ob- 
served ones if the spectra are rebinned to a res- 
olution of 70.000 (comparable to FTS spectra). 
Lacking visual phases for the models, the syn- 



' L.=7-103Lo, M. = 1.5Mo, r.=2600K, C/0=1.4, 
P=490d, AMp=6km-s"' 



thetic RVs in Fig.n(right) are shifted by an ar- 
bitrary value of 0.3. It can be recognised, that 
the results of the computations resemble the 
observations reasonably well (line-doubling, 
S-shape, amplitude ARV, etc.). This fundamen- 
tal characteristic of Miras can now be repro- 
duced even quantitatively in a self-consistent 
way by the new model. 
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